% Test the functions makeppftfilter and applyppftfilter against the
% preconditioned PPFT.
%
% Yoel Shkolnisky, December 2010.

for n=2:2:16;
    A=rand(n,n,n);
    tic;
    B=preconditionedAdjointPPFT3D(optimizedPPFT3D(A));
    t1=toc;
    
    H=makePPFTFilter3D(n,1);
    tic;
    C=applyPPFTFilter3D(A,H);
    t2=toc;

    err=norm(B(:)-C(:))./norm(B(:));    
    reportTestResult(n,err,t1/t2);
end
